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ABSTRACT 

Because of the corotation, the polarization angle (PA) curve of a pulsar lags the in- 
tensity profile by Ar / R\ c rad in pulse phase. I present a simple and short derivation of 
this delay-radius relation to show that it is not caused by the aberration (understood 
as the normal beaming effect) but purely by contribution of corotation to the electron 
acceleration in the observer's frame. Available altitude-dependent formulae for the PA 
curve are expressed through observables and emission altitude to make them imme- 
diately ready to use in radio data modelling. The analytical approximations for the 
altitude-dependent PA curve are compared with exact numerical results to show how 
they perform at large emission altitudes. I also discuss several possible explanations 
for the opposite-than-normal shift of PA curve, exhibited by the pedestal emission of 
B1929+10 and B0950+08. 
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1 INTRODUCTION 

In the simplest model of pulsar polarization (Komesaroff 
1970; Radhakrishnan & Cooke 1969, hereafter RC69) the 
position angle of polarization does not depend on the radial 
distanc^B of radio emission r. The polarization angle (PA) 
becomes dependent on r when dynamic effects of pulsar's ro- 
tation are taken into account (dragging of electrons by the 
corotating magnetic field). Blaskiewicz et al. (1991) (here- 
after BCW) showed that if the emission originates from a 
fixed radial distance r, the shape of the PA swing is (approx- 
imately) preserved (ie. it is the same as in the case of negli- 
gible r), but the entire swing is shifted towards later phases 
by A0bcw ~ <ir/R\ c radians with respect to the center of 
the profile (where Ri c is the light cylinder radius). Follow- 
ing BCW, I will refer to this formula with the term 'delay- 
radius' relation. Hibschmann & Arons (2001) (hereafer HA) 
have shown that the PA curve also undergoes vertical shifts, 
ie. in the PA values. Both these results appear to have in- 
teresting observational consequences (eg. Ramachandran & 
Kramer 2003; von Hoensbroech & Xilouris 1997). 

In the case of phase-dependent emission altitude differ- 
ent parts of the PA curve undergo different shifts and the 
PA curve assumes a distorted shape. This effect regularly 
happens to be employed to model observed distortions of 
PA curves and to derive magnetospheric emission altitudes 
(eg. Krishnamohan & Downs 1983, hereafter KD83; Xu, 
Qiao & Han 1997; Gil & Krawczyk 1997; Mitra & Seiradakis 



2004). A tool that is needed for this is an analytical formula 
for the PA that explicitly depends on the radial distance of 
the emission region. BCW and HA provide various forms of 
this equation. However, their formulae are not in a ready- 
to-use form: they are expressed through the emission time 
instead of the pulse longitude cf> b s (hereafter called pulse 
p/iasc0) . It is the need for this last step of the BCWs anal- 
ysis that actually sparked writing of this paper. 

A strict and formal description of the PA subject is 
given in the superb work of Hibschman & Arons (2001; see 
their appendices) and it will not be repeated in this pa- 
per. My intention here is to provide a simple reference for 
those who want to use the altitude-dependent PA curves 
in their data modelling. The aim is to clarify some obscure 
aspects of the subject by trivialising the formalism and to 
provide practical PA equations in their final form. Accord- 
ingly, Sect. [2] presents a very simple and short derivation 
of the delay-radius relation to clearly expose its origin. In 
Sect. |3] I introduce the fiducial phase, describe the magni- 
tudes and directions of various relativistic shifts with re- 
spect to it, and I write down the equations for the altitude- 
dependent PA curve in a form that is ready for immediate 
use in data modelling. In Sect. [4] I compare the analyti- 
cal approximations to exact numerical results obtained for 
various emission altitudes to show the validity range of the 
BCW theory. In Sects. [4]4]and[43]I discuss possible explana- 
tions for the opposite-than-expected shifts of PA curve, as 



1 The quantity r represents the distance measured from the cen- 
ter of the neutron star. 



2 Throughout this paper the phase is assumed to be measured in 
radians whenever dimensionless terms are added to it. 
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exemplified by the pedestal radio emission components of 
B 1929+ 10 and B0950+08. These and other interpretations 
of the anti-BCW shifts are summarized in Sect. [5] 



2 SIMPLE DERIVATION OF THE 

DELAY-RADIUS RELATION DUE TO 
COROTATION 

The derivation presented here is constrained to the case 
of the equatorial 3 ] plane of an orthogonally rotating pulsar 
(with dipole inclination a — 90°). However, it is simple, in- 
tuitive and demonstrates the effect more directly than the 
original derivation of BCW. 

Because of their corotation, trajectories of electrons are 
bent forward (toward the direction of rotation) in the in- 
ertial observer frame (IOF) with respect to trajectories in 
the corotating frame (CF). The trajectory of electrons that 
move along the dipole axis acquires some forward curva- 
ture, whereas the rectilinear motion occurs somewhere on 
the trailing side of the dipole axis. The bundle of IOF- 
trajectories becomes approximately symmetric around the 
location of this straight trajectory, which now assumes the 
role that the dipole axis had in the RC69 model (RVM 
model). This zero-curvature trajectory is associated with the 
inflection point of the PA curve. Its location can be found 
in the following two steps. First, we calculate 'rotationally- 
induced curvature', which is understood as the curvature 
imposed by rotation on a trajectory that is a straight line 
in the CF. Second, we will search for a place on the trailing 
side of the dipole axis, where this rotational curvature is ex- 
actly cancelled by the same-magnitude backward curvature 
of dipolar magnetic field lines. 

Relativistic electrons (with the speed of essentially con- 
stant magnitude v ~ c) that move radially along the dipole 
axis in the frame rotating with the neutron star's angu- 
lar velocity Q, undergo acceleration a ~ 2Qc in the IOF. 
This result can be obtained in several elementary ways (see 
the exercise in Appendix [XJ . The particle trajectory along 
the dipole axis thus acquires radius of curvature p r t = c 2 /a 
which is: 

RlC /,N 

Ptt ^ — • (1) 

The index 'rt' is to remind that this is the curvature induced 
by rotation on electron trajectory that is radial (or almost 
radial) in the corotating frame. Near the dipole axis, the 
curvature radius of dipolar field lines is given by 

4 r 
3 sin a 

where r = \r\ is the radial distance of an emission point and 
9 is the angle between r and the dipole axis. 

We look for regions in the trailing part of the polar cap 
tube, where the rotationally-induced curvature of eq. (fTJ is 
cancelled by the curvature of B-field lines. By equating p rt 
with pb one can find an equation for the locations of straight 
portions of the trajectories: 



8 r 



(3) 



where the index 'zc' stands for 'zero curvature'. We can see 
that the locations of the zero-curvature depend on radial 
distance: 9 ZC cx r. 

In the dipole geometry, magnetic field lines at points 
with coordinate 9 are directed at angle 9t — (3/2)9 with 
respect to the magnetic axis. Therefore, the phase delay of 
emission from the zero-curvature regions with respect to the 
(same altitude) emission along dipole axis is equal to 

s 3 „ . r 



A0bcw — 9k 



Ru 



(4) 



Thus, the locally-straight portions of electron trajectory lag 
the dipole axis spatially by 9 ZC ~ (8/3)r/Ri c , whereas the 
tangent-to-B directions of radio waves emitted from these 
regions lag the direction of waves emitted along the dipole 
axis by the angle of 4r/Ri c radians. Note that until the very 
end of the derivation there was no need to explicitly refer to 
the effects of aberration and propagation time delays (APT) 
that noticeably increase complexity of more general analysis 
(see Appendix [Bl for the definition of the APT effects). This 
is because it is the rotational straightening of electron trajec- 
tories that is the essence of the effect. As long as one is inter- 
ested in emission from a fixed altitude, both the aberration 
and propagation time effects can be ignored, because the 
waves emitted along the dipole axis as well as those emitted 
from the zero-curvature regions are advanced in (absolute) 
phase by roughly the same magnitude (r/R\ c by aberration 
and another r/R\ c by propagation time). It is so because in 
the small angle approximation all the open field lines are 
basically orthogonal to the corotation velocity (for details 
see Dyks, Rudak & Harding 2004b, hereafter DRH). The 
APT effects only have to be included when there are alti- 
tude differences, or when one wants to know the magnitude 
of absolute delays with respect to the "fiducial phase" . This 
subject will be discussed in detail in the next section. 

According to eq. ([3]), the last open field lines, located 
near 9\ ~ (r / Ric) 1 ^ 2 , have the zero-curvature in IOF at 
radial distance 



q 

4 -0.14. 



(5) 



3 Unless specified otherwise, 'equatorial' refers to the rotational 
equator. 



r 

rTc 

Note, however, that at this altitude the inaccuracy of our 
approximation (linear in r/R\ c ) is not negligible (it is of the 
order of (r/R lc ) 1/2 ~ 40%). 

We can also determine the rotation period P — 2it/Q, 
for which the zero curvature occurs at the trailing part of 
the rim of the polar cap. Since 9 zc /9i ~ (8/3) (r /Pic) 1 / 2 = 
0.0386(r 6 /P) 1/2 , with r 6 = r/(10 6 cm) and P in seconds, 
one obtains: 

P~ 1.49- 10~ 3 s Re, (6) 

where Re = Pns/(10 6 cm) is the neutron star radius. Thus, 
it happens just in the fastest known millisecond pulsars 
(J17148-2446 with P = 1.396 ms, Hessels et al. 2006, or 
B1937+21 with P = 1.558 ms, Backer et al. 1982) that the 
'dipole axis' is shifted full way to the trailing rim of the polar 
cap (provided the pulsars have a — 90° , as their interpulses 
suggest). The effect is illustrated in Fig. 1 of Dyks & Rudak 
(2002), where a numerically calculated trajectory of photons 
in the CF (for P = 1.5 ms and a — 90°) neatly coincides 
with the last open field line that emerges from the trailing 
side of the polar cap. Another consequence is that the ef- 
ficiency of magnetic pair production in these objects is the 
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weakest near the trailing side of the polar cap, not at the 
polar cap center (see Fig. 7 in the last-mentioned paper). 



3 POLARIZATION ANGLE CURVE 

In the absence of the rotational effects the PA as a function 
of phase is given by: 



tjj — tan 1 



sin a sin (<p obs 



cos(<^ ob s — 4>o) cos £ sin a — cos a sin £ 



+^M7) 



where (j>ob B is the pulse phase with zero point defined arbi- 
trarily by the observer, <f>o is the pulse phase at which the 
line of sight lies within the (SI, ft) plane, the constant ipn is 
the position angle of the projection of the pulsar rotation 
axis on the plane of the sky, a is the angle between the 
angular velocity of pulsar rotation SI and the dipole mag- 
netic moment fl, whereas ( is the angle between SI and the 
unit vector of the line of sight h h s ■ The sign of the arctan 
term is for ip measured in the observers' conventional way, 
i.e counterclockwise in the sky (Everett & Weisberg 2001). In 
this simple geometrical model, the PA is determined purely 
by projection of magnetic field direction on the sky's plane. 
Therefore, the azimuth angle (p^ of the dipole axis (mea- 
sured from the plane which contains an observer and SI, 
see FigfTJ and the pulse phase </> b s are simply related by 
<^m ~ ^obs — 4>o- Let us denote the value of pulse phase at 
which the center of the PA curve (its inflection point) is ob- 
served by 0pa (at this phase \dip / dcj>obs\ has a maximum and 
ip ~ V'n)- In the absence of the rotational effects, (and in the 
case of infinite propagation speed c = oo), the center of the 
PA curve would be observed at the pulse phase 4>pa = 4>o- 
Let us use the symbol <f> pIO { to denote the center of the pulse 
profile defined as the midpoint between the profile's outer 
edges. As long as the outer boundary of the radio emission 
region in the CF is symmetric with respect to the (fi, ft) 
plane, the center of the pulse profile would be observed at 
exactly the same phase pro f = 4>o- 

3.1 Altitude-dependent polarization angle curve 

BCW generalized eq. Q to include rotational effects. By 
considering the radio emission from a fixed radial distance r, 
they found that the PA curve is shifted toward later phases 
by A(/>bcw ~ ir/Ric with respect to the center of the pulse 
profile. For this special case of fixed r, the profile's center is a 
very convenient reference point from the observational point 
of view: if measured, the shift can be directly translated into 
the radial distance r/R\ c « A^bcw/40 For other applica- 
tions of BCW theory, however, the center of the pulse profile 
is not a suitable reference point, because the phase at which 
it occurs (0prof) depends on r itself. For example, the shift 
(PA center - profile center) does not tell us by how much 
the PA curve is shifted with respect to the case with the 
rotational effects ignored. Neither is the profile center help- 
ful if different parts of the profile originate from different 
altitudes. 



4 To find the phase of inflection point 4>VA one can fit the stan- 
dard PA swing (cq.[7} to the data (as BCW and von Hoensbroech 
& Xilouris did), because it has nearly the same shape as the equa- 
tion which includes the rotational effects. 




Figure 1. Orientation of the CF emission direction k 1 , IOF emis- 
sion direction k, and the dipole axis fl at the moment when ra- 
diation emitted at angle p with respect to fi in the CF becomes 
directed toward an observer who is located at angle £ from the 
rotation axis in the (Sl,x) plane. The line of sight direction n t, s 
coincides with k. Up to the order of r/R\ c the angles £ and Q 1 can 
be considered equal. The azimuths </>cm and </>' em differ by r/R\ c . 
Note that the angles (j> em and <j>' Bnl are assumed to increase in the 
opposite direction than </> M . 



A good reference point for measuring altitude- 
dependent shifts in pulse profiles should be unambiguously 
associated with the azimuth M of the magnetic dipole axis. 
This criterion is met by the fiducial phase 0/, defined as 
follows: it is the pulse phase at which the observer detects 
a photon that was emitted from the center of the star at 
the moment when the dipole axis was in the plane contain- 
ing fi and the observer (hereafter (SI, n b s )-plane). Because 
0obs — 4>f corresponds to M = strictly in the just- described 
sense, the fiducial phase <f)f plays a similar role as 4>o does 
in eq. [7](ie. one may consider them identical: <pf = <j>o). 

The shape of emission region that I consider in this pa- 
per is assumed to be axially symmetric around the dipole 
axis, ie. its radial distance r can be described as a function 
of only the angle p fa 1.5s(r/7?i c ) 1//2 between the emission 
direction in the CF and the dipole axis (s = SUT f/8 pc is 
the footprint parameter of a B-field line on the star sur- 
face). Such an average form is roughly consistent with the 
observed shapes of phase-averaged pulse profiles (eg. John- 
ston et al. 2008; Rankin 1983) but neglects the azimuthal 
structure (Karastergiou & Johnston 2007; Rankin & Ra- 
machandran 2003). To calculate the PA curve for an arbi- 
trary shape of r(p) one needs an equation for the altitude- 
dependent PA curve. Equations (16) or (17), (the latter with 



© 2002 RAS, MNRAS 000,[TJfT2] 



4 J. Dyks 



<f>o — 3r/R\ c ), from BCW can be used for this purpose, pro- 
vided that a transition from 'fit' to the pulse phase <p h B is 
carefully done. 

The equations (16) and (17) of BCW are expressed in 
the emission tim«0 with the zero point defined in a particu- 
lar way: t = t e = corresponds to the moment when the 
dipole axis is in the (fi, n b a )-plane. Thus, the emission time 
in BCW is simply defined as t e = 0^/fi, where 4>n is the 
azimuth of fl in the frame with fi \\ z and with the observer 
in the (fi, x) plane (see Fig. [T] this azimuth provides the 
measure of the emission time). It is very important to dis- 
cern the emission time (or </> M ) from the detection time td 
(or from pulse phase (pohs) at which the radiation emitted 
at t e is detected. They are related by: 



td = t e + d/c - r ■ n obs /c -)- At zp , 



(8) 



where n b s is the unit vector pointing towards the observer, 
d is the pulsar's distance and At zp takes into account the fact 
that the observer is allowed to assume arbitrary zero point 
in the counting of time. The term f • h \, B /c ('propagation 
time advance') takes into account the fact that the source 
located at f is closer to the observer than the center of the 
neutron star. The same can be expressed in terms of angles: 



(f)^ + fid/c — fir ■ hobs/c + Ac6 zp , 



(9) 



where ci bs = fitd, and A<j> zp takes into account the fact that 
the observer can assign phase zero to the pulse profile in 
an arbitrary way. The pulsar distance d and the zero point 
difference Acp 2p can be replaced in this equation with the 
fiducial phase. The definition of <pf, as articulated above, is 

<f>f = 0obs(0 M = °, r = °) = Orf/c + A^ Z p (10) 

(from eq. [9}, which can be inserted back into (0 to obtain: 

<^obs =4>p — Qr- n obs /c + <f>f w <£ M - r/Ri c + (f> f , (11) 

where on the right hand side I make the small-angle approxi- 
mation: f -n b s ~ r and use fi/c = 1/R\ C . Remembering that 
<ji M = fit e we get: 



r/Rl, 



(12) 



which should be used in equations (16) and (17) of BCW 
in place of 'fit' (the common practice of replacing fit in 
eqs. (16) and (17) of BCW with <^ b s results in PA curves 
that lag the fiducial phase by 3r/Ri c instead of the actual 
2r/R\ c ; therefore they underestimate r by a factor of 1.5.) 

Thus, the equation for the altitude-dependent PA curve 
reads: 



Ipr 



tan 



3(r/Ri c ) sin £ — sin a sin( 



■<p f +r/R u 



sin P + sin a cos ((1 — cos[</> b s 



<pf + r/R lc ]) 
(13) 



+ 



where /3 = £ — a is the observer's 'impact' angle. The equa- 
tion was obtained directly from eq. (16) in BCW (in addi- 
tion to the use of fit = <ji> bs — <pl + r/Ri c I have changed the 
BCWs sign of the arctan term to agree with the observers' 
convention) . 



5 The emission time refers to a Lorentz frame in which the pul- 
sar's center of mass is at rest. Note that BCW use two different 
symbols to denote it (t as well as t e ). 



The same approach may be applied to eq. (17) of BCW 
(with 4>o — 3r/Ri c ) to get: 



! tan 



- sin q sin^0 obs 



2r/R h 



sin/3 + [1 - cos(</> obs 



10 r 

+ ^--5-coscv-|-V«, 

6 Hie 



2r/Ri c )] cos £ sin a 
(14) 



which is equivalent to (|13[) within the accuracy of the BCW 
method, ie. up to the order of r/R\ c . The manually-added 
term (10/3)(r/7?i c ) cos(a) represents the vertical shift of the 
PA found by HA (the shift given in HA refers to the clock- 
wise definition of PA, and, therefore, has opposite sign). For 
the polar-current flow that is close to the Goldreich- Julian 
value, this term is cancelled and can be neglected in eqs. (|14|l 
and (| 1 5f) . In such a case, however, the current-induced shift 
of PA must also be subtracted from eq. (|13[1 for consistency 
(see Sect. 14.21 below). 

Yet another, most direct method to derive the equa- 
tion for the altitude-dependent PA curve is the following. 
The center of the pulse profile precedes in phase the cen- 
ter of the PA swing by 4r/R\ c , ie. 0pa ~ pro f + 4r/Ri c . 
The profile center itself precedes the fiducial phase <j>f by 
2r/R\ c (one r/R\ c for the aberration and another r/R\ c for 
the propagation time, see eg. DRH for details). Therefore, 
the center of the PA curve lags the fiducial phase </>/ by 
Ar/Ri c — 2r/R\ c — 2r/R\ c and at any phase 4> hs the posi- 
tion angle ip r which includes the rotational effects is equal to 
the 'nonrelativistic' PA (given by eq. taken at the earlier 
phase </>obs - 2r/R\ c . Thus, Vv(</>obs) = ip(<f>ob a -2r/Ri c ): 



Ipr 



' tan 



sinasin(c/> b s - (pf - 2r/R lc ) 



COS(0 obs 



10 r 

+ — ~fz — cos a + ipn 

6 Kic 



2r /Ric ) cos £ sin a — cos a sin £ 
(15) 



This equation was simply obtained by replacing o b s in 
eq. with o bs — 2r/R\ c and by replacing <po with <pf. 
Simple trigonometry shows that it is equivalent to eq. (|14[) . 

In the limit of r < 10" 2 i?i c all eqs. (JT3J - (JT5j) reduce 
to eq. (0 and the centers of PA the curve, the pulse pro- 
file, and the fiducial phase coincide. In eqs. (I13p - (|15|l the 
radiation electric field E w is assumed to be along, rather 
than orthogonal to, the direction of electron acceleration a 
(in eq.[7]i5 u , is assumed to be along the magnetic field B, or 
at least in the plane of a £?-field line). For the case E w _L a 
(or for E w LB), the value of PA given by eqs. (0, (fT3)l . 
(|14|) . and (| 15|) may need to be increased by 90° (as in the 
case of the Vela pulsar, Lai et al. 2001, Radhakrishnan & 
Deshpande 2001). Obviously, if the projection of the rota- 
tion axis at the plane of the sky is unknown the terms ipa in 
eq. (0 and (|13[) - (|15[) must be considered free parameters 
to be determined from the fit to the data. If r is independent 
of pulse phase and we are not interested in the value of ipn, 
the term (10/3)(r/-Ri c ) cos a can be merged with ipn into a 
single fitting parameter. 

Eq. (|15p depends only on the term —<pj — 2r/R\ c which 
makes it difficult to separate <j)f from r. A special case when 
it is possible is when r does not change across the pulse pro- 
file. In such a case the center of the PA curve (the inflection 
point) is located at the phase 



6pa ~ 4>f + 2r/R\ c . 



(16) 
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The center of the pulse profile is then located at 

<?W ~<Pf- 2r/R lc (17) 

so that the shift of the PA with respect to the profile is 
4>pa — 4> P rof ~ 4r/7?i c . Then (i.e. for r — const) eqs. (|16[) 
and |17J) tell us that 

0/ = (&>rof +4>pa)/% (18) 

ie. 0/ is half way between the center of the pulse profile and 
the center of the PA curve (Fig. [2]). If the nature and the 
quality of the data allow us to determine </> pro f and </>pa , the 
value of 4>f can be calculated from (|18|) under the asumption 
that r = const. In this specific case <^>pa can be determined 
by fitting the nonrelativistic PA curve (eq.O to the observed 
data. 

3.2 Relation between the radial distance of radio 
emission and the pulse phase at which it is 
detected 

In the case of emission altitude that changes gradually with 
pulse phase </> bs, the eqs. (|13[1 - (|15|l must be supplemented 
by an analytical equation for r — r(4> b B ) to be useful. 

If the distribution of emissivity in the CF frame is sym- 
metrical with respect to the dipole axis, the radial distance 
of the radio emission r is a function of only the angle p be- 
tween the emission direction in the CF and the dipole axis. 
Let us assume that the emission region can be described by 
a simple function r(p) which can be invertecjf] into p(r). 

Fig. [T] shows relative orientations of the dipole axis p, 
and the CF-emission direction k' at the moment when the 
radiation is directed towards the observer. In the IOF it 
propagates along k \\ n b B - 

The radiation from some point at radial distance r is 
directed towards the observer only when p is rotated by an 
appropriate angle <^) M = <f) CYn — Qt em . Because the aberration 
advances the radiation by r/Ri c in phase, the angle cm is 
smaller by r/Ri c than the angle by which the dipole 
would need to be rotated in the absence of the aberration 
to become aligned with the line of sight (see Fig. [TJ: 

4>' cm ~ em + r/Ri c = tttcm + r/R lc . (19) 

From (|12[) and (|19[) we learn that the radiation is shifted 
toward earlier phases by 2r/R\ c with respect to <f>j and is 
detected at the phase: 

obs ~0 / + 0e m -2r/iii c . (20) 

By substituting (j>ob B in eq. (|15[1 by the above formula, and 
using 4>' crn ~ cm + r / Ri c , as well as <f> cm = Q,t cm one can 
easily verify that eq. (I15|l is equivalent with eq. (17) in BCW. 

Unlike (j> eul (see Fig. [TJ, the azimuth (j>'cm of the non- 
aberrated emission direction k' is associated with the angle 
p between p and k' through: 

cos [p(r)] = cos 4>' eul sin a sin £ + cos a cos f (21) 




-15 -10 -5 5 10 15 

pulse phase [deg] 



Figure 2. Influence of increased emission altitude on the observed 
PA curve. Panel a) shows the RC69 case of negligible emission 
altitude. In b) the emission altitude has been rised up to 0.03i?i c . 
As a result, both the intensity profile and the PA curve move 
in opposite directions from the fiducial phase by the same angle 
of 2r/R\ c rad. In both panels pulse phase zero is the fiducial 
phase ('dipole axis phase'). In BCW, the only figure that has the 
fiducial phase at </> b s = is their fig. 3. The apparently small bar 
at o bs = 10° shows the vertical upward shift of the PA curve 
by (10/3)(r/iJi c ) cos a. The solid PA curve in b) presents the 
approximate equation of BCW (eq. I13l in this paper). The dashed 
line is just the appropriately shifted RVM curve (eq. I14jl , In a) 
they coincide. In the figure a = 45° and £ = 41.3°. In all figures 
V>n =0. 

which is a direct form of the spherical trigonometry cosine 
theorem applied for the triangle (Q,k',p). Calculating <f>' em 
from (|21|) and inserting into eq. (|20[) gives 

, x -if cos(p(r)) — cos a cos C\ , „ r ,„„, 

0obs r ~±cos — i +^-2—, (22) 

V sm a sin Q J Ri c 

where we ignored the insignificant difference between £ and 
The '+' sign at the arccos term corresponds to the trailing 
whereas the ' — ' sign to the leading part of the open field line 
region. 

Since the argument of the function arccos in eq. (|22|l 
does not have to be small (e.g. for small dipole inclination 
a the profile width can reach several tens of degrees), the 
equation in general cannot be inverted to obtain a simple 
analytical formula for r(<j!> bs). To use eq. (|15[1 for phase- 
dependent r, the latter needs to be determined from (|22|l 
numerically. 



6 If more than one value of r correspond to the same value of 
p (ie. if there are several layers of emission located above each 
other) one can separate r(p) into a few functions p, each of which 
is reversible. Since emission from different layers can be observed 
simultaneously, the Stokes parameters must be used as eg. in 
KD83 or Mitra & Seiradakis (2004). 



3.3 Numerical example 

Eq. (|22[) is useful for a quick examination of the shape of 
PA curves for a variety of emission regions, ie. for various 
functions r(p). A convenient way of doing this is to define 
a dense table of angles pi, and calculate the corresponding 
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— = 0.01 + 0.02 sin , (23) 

Ric \ ao ) 

where the constant ao = 2° was used in panels a and b, 
whereas ao = 1° in c and d (this somewhat strange shape 
was chosen to make the distortions of the PA curves easily 
noticeable by eye). The curves r(0 o b s ) (as given by eq. [22)) 
are shown in panels b and d. One can see the characteristic 
skewing of the r(0 o b s ) function toward early phases, which 
is caused by the APT effects. In panel d, the radial dis- 
tance r increases with phase so fast that the radiation from 
the upper parts of the emission region overlaps in phase 
with radiation emitted at lower altitudes. The correspond- 
ing PA curve (panel c) is not a mathematical function of 
phase (there are a few regions with three values of PA re- 
ferring to the same o bs). To look normal it needs to be 
cut into pieces and Stokes-summed (for details see Mitra & 
Seiradakis 2004; KD83). 



15 



: c) 




— i — i — i — i — 

id) \ 


I — i — i — i — i — h 


— i — i — i — i — | — i — r 


— i — i — | — i — t 


— i — i — I — ' 





15 



Figure 3. Altitude-dependent PA curves (thick lines in a and 
c) calculated by using eqs. i ]22[ ) and (1 1 3 I t for the emission region 
of eq. J23} with a = 2° (top) and a = 1° (bottom). The thin 
line shows the standard RVM curve for reference. Panels b and 
d show the corresponding radial distances in units of R\ c . Phase 
<t>obs = is the fiducial phase and I used a = 45° and £ = 43° in 
the figure. 



vector of rt(pi) for some chosen function r(p) (or vice versa: 
to define rt and calculate pi{ri) for arbitrarily selected p(r)). 
These are next used in eq. (|22p to calculate the table of the 
corresponding values of o bs,i- The tables rt, and </> b s ,i can 
then be directly used in eq. (|13p (or I25p to calculate the 
shape of the PA curve ip(4>obs,r). 

Fig.|3]presents an example of such procedure performed 
for the emission region given by: 



4 LIMITATIONS ON THE VALIDITY AND 

APPLICABILITY OF THE LINEAR THEORY 

The equations of Sect. [3] refer to phase shifts that are usu- 
ally tiny and can be easily affected by several interfering 
effects discussed in this section. Moreover, the accuracy of 
the analytical results of previous sections is limited by their 
inherent error of (r/Ric) 1 ' 2 . The validity of the analytical 
theory is therefore severly limited: on one hand by the too 
small magnitudes of the phase shifts (too small to be re- 
liably measured, or to be unaffected by the disturbing ef- 
fects); on the other hand by the poor accuracy of the theory 
at larger r. The following subsections are to show that it 
makes practically no sense to apply the theory for r > 0.1J?i c 
(A0 obs > 10° -20°). 



4.1 Asymmetry of the open field line region 
caused by rotational distortions of the 
magnetic dipole (sweepback) 

The strict symmetry of the open volume with respect to 
the (Q, /2)-plane is unlikely, and it has been shown (Dyks & 
Harding 2004, hereafter DH04) that even very small effects 
(that at low altitudes are of high order in r/R\ c ) can eas- 
ily introduce large asymmetry. Here I discuss the rotational 
deformations of the vacuum dipole (DH04; Shitov 1983, 
Deutsch 1955), because they are easy to analyze (thanks to 
the existing analytical formulae) and they seem to be qual- 
itatively similar to those obtained in recent plasma-loaded 
simulations (eg. Spitkovsky 2008). As we show in detail in 
DH04, in the lowest order (of (r/R\ c ) 2 , which is already 
higher than the r / R\ c accuracy of equations in the preced- 
ing sections) the rotational sweepback does not introduce 
any asymmetry to the B-field direction. The asymmetric 
change of B-field direction is of the order of (r/i?i c ) 3 at 



low altitudes (completely negligible there). Despite this, a 
very strong asymmetry of the open volume (of magnitude 
~ (r / Ric) 1 ^ 2 ) is generated by the sweepback. This is be- 
cause the low-altitude outer boundary of the open volume 
is determined by the geometry of B at the light cylinder, 
where the sweepback becomes a strong effect (first order, in 
the sense that there we have (r/R\ c ) a ~ 1, regardless of a). 
This strong distortion is traced back towards low altitudes 
through the continuity of magnetic field lines, and results 
in strong (much stronger than (r/Ri c ) 3 ) asymmetry of the 
open volume. 

This effect can be taken into account by replacing the 
delay-radius relation (eq.|3| with the 'misalignment' formula 
of DH04: 



A0~4-^- -f( 
v Ric \R lc J 



1/2 



(24) 



where F depends on a and C, and has typically the magni- 
tude of ~ 0.2. Eq. (|24[) differs from the original delay-radius 
relation only in that it takes into account the azimuthal 
asymmetry of the outer boundary of the open field line re- 
gion. However, I do not refer to eq. (|24[) with the name 
'delay-radius' relation because the shift of position angle 
curve does not have to be the positive delay anymore - it 
can now slightly precede the profile midpoint, see DH04 for 
more details. 

Those, who believe that the rotating vacuum dipole is a 
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Figure 4. Variations of polarization angle curve due to increasing emission altitude, marked in the bottom right corners of each panel. 
The dots, often merging into a thick solid line, present exact numerical solution. The thin solid line presents the equation H13I I. The 
dashed line is for eqs. (114 I t and/or dl5l . For each altitude the numerical curve is shown for both the static-shape dipole (letter-marked 
panels) and the rotationally-distorted vacuum dipole (below the static case). The figure was calculated for a = 45°, £ = 41.3°, 4>f = 0, 
and ip n = 0. 



better approximation of the real B than the static vacuum 
dipole should use eq. (|24|l instead of the original delay-radius 
relation of BCW. What may be most reasonable is to use 
both the eq. (|24l) and ((4j as the lower and upper limit for r, 
respectively. 



4.2 Current-induced distortions of magnetic field 

The magnitude of the current-induced distortions of mag- 
netic field can be easily estimated with a simplification of 
Ampere's law V x B — c _1 (47rj + dE/dt) into the crude 
form of AB/L ~ J/c, where AS is the magnetic field gen- 
erated by the current density J, and L is the characteris- 
tic scale. For longitudinal polar currents of the Goldreich- 
Julian (GJ) magnitude we have J ~ Jgj ~ Pgjc, where 
pGJ — Q.B/(2nc) ~ B/Ric is the GJ density, whereas the 
transverse scale L ~ r9\ a ~ r(r/i?i c ) 1//2 corresponds to the 
open field line region. The distortions of B are therefore of 



the order of (AB/B) po i ar ~ (r/R lc ) 3/2 . For toroidal cur- 
rents (due to the corotation of Pgj) we have J ~ poj^rot 
where i> rot = Qr = cr / R\ c is the local corotation ve- 
locity, and L ~ r. This results in smaller distortions of 
(AB/B) tldl ~ (r/R lc f. 

HA find that the polar currents shift the PA curve 
downward by Aipj = — (W/3)(r / R\ C )(J/ Jgj) cos a. A spec- 
tacular way of estimating the current density has recently 
become possible thanks to the observations of part-time pul- 
sars (Kramer et al. 2006) and suggests that J ~ Jgj during 
the 'on' phase of pulsar emission. 

The simplest way to include the effects of current in 
the approximate eqs. (|13|l - (|15|l is to subtract the term 
(10/3)(r/_Ri c )(J/JGj) coscv with J/ Jgj = 1 from their right- 
hand sides, to obtain: 



1p r 



tan 



3(r/R\ c )smC - sinasin(0 obs — <j>j + r/R\ c ) 



sin/3-|- sinacos((l — cos[0 obs - 4>f + r/R\ c ]) 



+ 
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10 r / J \ 

for the BCW formula of eq. (11^1 , and 

sinasin(0obs — <j>f — 2r/R\ c ) 



(25) 



i tan 



cos(</) bs — 0/ — 2r/R\ c ) cos £ sin a — cos a sin £ 
+ i>n (26) 
for the 'shifted RVM' curve of eqs. (JT4J) and (fT5)l . 

4.3 Shapiro delay 

To make things possibly simple, the importance of general 
relativistic time delays has been neglected so far. These 
should not affect the delay-radius relation, because it refers 
only to a single r, and is insensitive to any time delays 
(whether gravitational or flat spacetime). They do affect, 
however, all formulae which refer to altitude differences (as 
do all the equations that explicitly involve both r and </>/). 
Let us then use the Schwarzschild metric to learn when the 
Shapiro delay cannot be neglected. The time for a ray to 
travel radially in the gravitational field from SOII16 f*min lip 

to fniax IS 



Ats 



dr 



c(l-r g /r)' 



(27) 



where r g — 2GMns /<? . Instead of that, our flat spacetime 
formulae assume the time is the Romer delay (r max — r mm )/c, 
so they miss the fact that the low altitude emission compo- 
nents (emitted at r m i n ) are additionally Shapiro-delayed by 
Ai Sh = Atschwrzid - (r ma x - r min )/c with respect to the 
high-altitude emission components emitted at r max . 

In the case of millisecond pulsars (MPs) this effect can 
become easily measurable: for r m i n = 10 6 cm, r max = 0.2i?i c , 
and P = 3 ms, we obtain Atsh = 2-10 -5 s which corresponds 
to the phase shift of 2.48°. For the same r m i n and r max /i?i c 
but for P = 1.5 ms the phase shift is 1.95° . It become slightly 
smaller because for the fixed r max /_Ri c = 0.2 the altitude dif- 
ference decreased (for P — 3 ms R\ c is 14.3 • 10 6 cm, whereas 
for P = 1.5 s it is twice smaller). Note that in the estimate 
I take r max = 0.2i?i c as a blind guess of the upper limit for 
radio emission regions in millisecond pulsars. This is simply 
because the known estimates or radio emission altitudes for 
classical pulsars (a few tens of Rns, DRH; Gupta & Gan- 
gadhara 2003; Kijak & Gil 1997) would locate the emission 
region beyond the light cylinder in the millisecond pulsars. 
If one takes r max = Ri c , the gravitational phase shifts be- 
come 5.44° and 8.41° for P — 3 and 1.5 ms, respectively. 
Thus, the equations that rely on altitude differences cannot 
be reliably applied to fast millisecond pulsars (in particular, 
it is not safe to use eq. (|17|l to determine emission altitudes 
from the core-cone shift). Another reason for not applying 
the approximate theory to the MPs is that it is very in- 
accurate for altitude differences of Ar > 0.1i?i c , which are 
likely in fast MPs even if the radio emission extends radi- 
ally by only one 7?ns- To calculate the phase shifts for MPs 
even more accurately, one would have to use the metric for 
a fast rotating neutron star (eg. Braje et al. 2000; Gonthier 
& Harding 1994; see also D'Angelo & Rafikov (2007) for 
a detailed description of how the altitude differences affect 
timing). 



+ 



In the case of classical pulsars (with P ~ 1 s and r max ~ 
50i?Ns), the Shapiro delays can be safely ignored, because 
they comprise a tiny fraction of the (now longer) rotation 
period. The phase shifts typically have a few hundredths of 
a degree. 



4.4 Limited accuracy of the lowest-order theory 

The relativistic phase shifts are of the order of r / R\ c , 
whereas the lowest-order terms that have been neglected 
in derivation of eqs. (|13|l - (|15p have magnitude of ~ 
(r/i? lc ) 3 / 2 . This means that the analytical approxima- 
tions for the polarization angle have a fractional error of 
100(r/i?i c ) 1//2 percent, which is as large as 30% already at 
0.1-Ric- The approximate equations (|13[) - (|15p are compared 
to the exact numerical PA curves in Fig. [4] The numeri- 
cal PA curves were calculated for fixed emission altitudes 
(shown in bottom right corners) and only involve the kine- 
matic effects of corotation (the method of calculation is de- 
scribed in sect. 2 of Dyks et al. 2004a). It is seen that for any 
r > 0.1-Ric there are considerable differences between the ex- 
act result (dots) and the approximate formulae. One can also 
see that at large altitudes the approximate formula of BCW 
(eq. 1131 thin solid lines) performs much better than the ap- 
propriately shifted classic equation of Komesaroff (eqs. 1141 
and 1151 dashed lines). 

An interesting effect seen in Fig. U is that for large al- 
titudes the numerical PA curves strongly tend to assume 
the distorted-sine-like ('equatorward') shape, in spite of that 
the figure is for the poleward viewing geometry (a = 45°, 
£ = 41.3°). The range of £ with the equatorward PA shape 
is not limited to the range (a, n — a) and is instead increas- 
ing with altitude. For the specific case shown in Fig. U the 
numerically determined range of the equatorward PA curve 
is (43°, tt - 43°) for r = O.LRic, and (38.5°, tt - 38.5°) for 
r — 0.2i?i c (with numerically-limited accuracy of ~ 0.5°). 
The BCW equation for the PA (eq. I13|l does reproduce this 
behaviour (to some limited degree), but the fixed-shape PA 
curve of the 'relativistically shifted' RC69 model (eq. Q3] or 
I15[) fails to do this, which results in the large disagreement 
already visible in panel c of Fig. [4] (for r — 0.2i?i c ). 

Another interesting effect in Fig. [4] can be seen if one 
compares PA curves of the static-shape dipole to the dis- 
torted ('retarded') dipole case. At large altitudes (r > 
0.7i?i c ) the sine-like PA curves for the distorted dipole 
(Fig. HJiji) move leftward (toward earlier phases) and have 
their steepest gradient point shifted leftward with respect to 
the fiducial phase (opposite than expected from the BCW 
theory). This effect is caused by backward bending of mag- 
netic field lines (the sweepback) which shifts the spatial lo- 
cation of "magnetic axis" (here I mean a locally straight 
B-field line at a given r) foreward, into the leading part of 
pulsar magnetosphere. At large r this effect seems to domi- 
nate the straightening of electron trajectories. 

The distorted-sine shape of the high-altitude PA curves, 
with the steepest gradient point well ahead of the 'dipole- 
axis' phase, bears close resemblance to the unusual PA curve 
of the pedestal radio emission of B1929+10, where the steep- 
est gradient precedes the main pulse by as much as 18° 
(Rankin & Rathnasree 1997; Everett & Weisberg 2001). This 
interpretation would imply that the pedestal radio emission 
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originates from r ~ 0.7-Ri c . However, below I also mention 
other possible interpretations of this effect. 

4.5 Limitations due to the specific radio emission 
mechanism assumed in the theory 

The radio emission mechanism assumed in the BCW theory 
is the curvature radiation caused by the 'macroscopic' ac- 
celeration due to the curvature of magnetic field lines. The 
acceleration is small enough that the corotation can modify 
it to produce the delay-radius relation of eq. @. The oper- 
ation of the curvature radiation in pulsar magnetosphere is, 
however, still a matter of debate (eg. Luo & Melrose 1992; 
Kunzl et al. 1998). 

The other emission processes (eg. the direct or inverse- 
Compton-scattered plasma emission, or the synchrotron 
emission) typically involve much larger accelerations due to 
microscopic motions, that are unlikely to be noticeably af- 
fected by the corotation (eg. the acceleration due to gyration 
exceeds the macroscopic acceleration due to the corotation 
by a factor cj B /tt ~ 10 8 [B/(10 6 G)](7/10 3 ) -1 , where uj b is 
the gyration frequency, Takata et al. 2007). 

It is then tempting to speculate that the PA curves 
for such emission processes would not be delayed by the 
BCW effect. But the aberration and propagation time effects 
would still be there, and would shift the PA curve by 2r/R\ c 
toward earlier phase than the main pulse. This may be a 
mechanism responsible for the leftward shift of the pedestal 
PA in B1929+10. In Dyks, Rudak & Rankin (2007; hereafter 
DRR) we provide other arguments for non-curvature origin 
of the pedestal emission. The main pulse emission, on the 
other hand, does not exhibit the abnormal shift which makes 
it more consistent with the curvature radiation. 



5 ANTI-BCW SHIFTS OF PA CURVE 

The best examples of this phenomenon are provided by the 
pedestal radio emission of B1929+10 and the weak bridge of 
radio emission that connects the main pulse and interpulse 
of B0950+08. The PA curves of these extended components 
(fitted without the part under the main pulse, see Everett 

6 Weisberg 2001) have their steepest gradient points well 
on the leading side of the main pulse. So far the following 
mechanisms have been considered for this effect: 

(i) The high-altitude curvature emission from roughly 
fixed altitude (Sect. 14.41 in this paper and Fig. HJi,i)- This 
probably cannot explain the case of B0950+08, for which 
the PA curve has the poleward shape. 

(ii) Radio emission mechanisms different than curvature 
radiation (sect. [43}. Possibly consistent with the parallel- 
ICS interpretation of double notches (DRR). 

(iii) Inward curvature radiation. Proposed simply as an 
inverted version of the BCW effect (Dyks et al. 2005). The 
'pulsar shadow' model of double notches, that has led to this 
interpretation, has been superceeded by the model described 
in DRR. 

(iv) Yet another interpretation can be devised using the 
'limiting polarization radius' idea: polarization of the ra- 
diation propagating through the inner (dense) regions of 
the pulsar magnetosphere can follow the local direction of 



magnetic field (eg. Cheng & Ruderman 1979; Melrose 1979; 
Barnard 1986) until the radiation reaches the polarization- 
limiting radius r po i. The plasma-density-dependent esti- 
mates of r po i are highly uncertain (eg. Lyubarsky 2002) but 
it should be significantly larger than the radial distance of 
the emission region. 

High-up in the pulsar magnetosphere the polar beam of 
radio emission, emitted at low altitudes, propagates through 
the trailing part of open field line region (see fig. 1 in Dyks 
& Rudak 2002, or fig. 1 in Lyubarsky 2002). Therefore, 
the polarization imprinted in the beam at r po i reflects the 
trailing part of the polarization curve, with the steepest 
gradient point shifted considerably towards early phases. 
The same mechanism was proposed by Lyubarsky (2002) 
to explain the shallowness of PA curves observed in mil- 
lisecond pulsars (Xilouris et al. 1998). Let me emphasize 
that the parallel-ICS model of pedestal emission mentioned 
in point (ii) above, actually necessitates this propagation- 
induced linear polarization. Otherwise the observed polar- 
ization degree would be very low due to the convolution of 
many micro-beams contributing to the line of sight. 



6 DISCUSSION - ABERRATION OR NOT? 

It is shown in Sect. [2] that the magnitude of the shift given 
by the delay-radius relation has nothing to do with the aber- 
ration (understood as the normal beaming effect). On the 
other hand, Hibschman & Arons mention a phase shift of 
3r/Ri c (their eq. 4) and claim that it is just the aberration 
("simple beaming") that is responsible for the remaining 
r/Ri c part of the total delay. Do we have a contradiction 
here? 

Not really, though the wording of HA can easily be mis- 
understood. Both in the BCW and HA, the delay-radius re- 
lation is derived as a difference of two moments of emission: 
1) the moment when the radiation bearing the steepest- 
gradient PA is directed to the observer (this happens at 
flt e = 3r/Ri c ) and 2) an earlier moment when the radiation 
in the middle of the profile (emitted in the (H, /I)-plane in 
the CF) becomes directed toward the observer. HA explic- 
itly notice that the latter does not take place at the moment 
t e = 0, when the dipole axis is in the f2-observer plane: in- 
stead, because of the aberration it happens at £lt e = —r/R\ c 
(see Fig. [1] and eq. [19] with <f>' crn = 0). The key point that 
is not mentioned in their description, however, is that the 
moment £lt e — 2>r / R\ c at which the steepest gradient radi- 
ation is directed toward the observer has already been also 
advanced in time by the same aberration angle of r / R\ c . 
Thus, the aberration, understood as the normal beaming ef- 
fect, works identically both at the center of the main pulse, as 
well as on its trailing side and contributes practically noth- 
ing to the shift of the PA curve with respect to the intensity 
profile. 

Another issue related to aberration is mostly nomen- 
clatural. In this paper I describe the origin of the delay- 
radius relation as the 'straightening of electron trajectories' 
when they are transformed from CF to IOF. Or, one can say 
the effect is due to the transformation of electron accelera- 
tion from the non-inertial CF to IOF. In the wording of HA 
the phase shift caused by these effects is described as 'aber- 
rational' too (see eg. their Appendix F and G) . An argument 
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that can (possibly) justify this is eq. (|A1[I (their eq. Fl). The 
addition of velocities (aberration) given in this equation de- 
termines the electron velocity in IOF, which is next differ- 
entiated to obtain IOF acceleration. Of course, what really 
matters is not the fact that we add the corotation velocity, 
but the fact that the added velocity is time-dependent (or 
that CF is non-inertial) . Since the phenomenon of aberra- 
tion by itself refers only to the velocity of a reference- frame, 
and not to its acceleration, the term 'aberration' misses the 
essence of the effect, which is the time- dependent nature of 
the corotation velocity in ()A1|I . 

A word of comment on papers that have assumed that 
the PA curve is shifted forward in phase, just as the intensity 
profile does (eg. KD83; Xu et al. 1997; note that the work 
of KD83 was published in the pre-BCW era). If one corrects 
the analysis of KD83 for the direction of the PA shift, the al- 
titude order of their four emission regions (see. their fig. 16) 
should be turned upside-down, ie. the region 4 should be at 
the top, region 3 below 4, etc. 

Finally a note related to the high-energy emission. Elec- 
trons with Lorentz factors 7 ~ 10 6 - 10 7 that move along 
the dipole axis trajectory with p r t = i?i c /2 and P ~ 0.1 s 
emit curvature radiation that extends up to 0.1 - 100 MeV, 
respectively. This energy is too low for pair production (see 
fig. 7 in Dyks & Rudak 2002). However, this radiation may 
contribute to the observed X-rays and gamma-ray pulse pro- 
files and spectra. In the numerical models of the high-energy 
emission from pulsars it is therefore necessary to calculate 
the curvature radiation using the radius of curvature of elec- 
tron trajectory in the inertial observer frame. 
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APPENDIX A: SIMPLE DERIVATION OF 
EQ. (1) 

Let us consider the geometry shown in Fig. IA1I The three 
positions A', P', and B' of a point at the light cylinder are 
separated by the same path increment s ~ R\ c 5 so that the 
triangles OA'P' and OP'B' are identical. As soon as we mark 
in this figure some arbitrary position, say A, of an electron 
that is moving radially in the corotating frame, the trajec- 
tory APB of the electron becomes uniquely determined. This 
is because: 1) the electron also moves with the speed of light 
c (as the point at the light cylinder does), so that we must 
have \AP\ — \PB\ = s and 2) the electron has the same 
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Figure Al. Trajectory APB of a relativistic electron that is mov- 
ing along the straight dipole axis in the corotating frame. Rota- 
tion axis is perpendicular to the page at O. The points A', P', 
and B' present locations of a point corotating at the light cylin- 
der. The diagram is used in the school-style derivation of cq. ||T)| 
(Appendix lAl , 



angular velocity 17 so that both AP and PB must subtend 
their respective angles S. 

The electron trajectory has a radius of curvature p It = 
s/e (for elementary reasons), where e — n — k and k is the 
angle between AP and PB. For the quadrangle OAPB we 
can write 28 + u + k + u' = 2tt which implies e = 28 + u) + 
to' — 7r. The sine theorem for the OPB triangle gives sin uo' = 
rsin(5)/s — r8/s ~ r/Ri c , i.e. ui' ~ arcsin(r/i?i c ) ~ r/R\ c . 
For the OAP triangle one has smu — rsin(<5)/s ~ r/Ri c 
(same as for lj') but this time u) is in the range tt/2 < ui < n 
(see Fig. IA1|) so that ui = ir — arcsin(r/i?i c ) ~ n — r/R\ c . 
Thus we have e ~ 28 ~ 2s/R\ c so that p rt = s/e ~ R\ c /2. 

A formal derivation of eq. Q assumes that the electron 
velocity is given by: 



cb + fi x 



(Al) 



where 6 is a unit vector along the local direction of B, and 
the tiny difference between the electron speed along b and c 
is neglected (hence the '~' sign). The electron acceleration 
is then: 



dv db =< dr dl > 

a — — = c— + il x — — c— + il x v 
dt dt dt at 

db 



(A2) 



c— +c(v- V)b + cfl x 6 + x (fi x r), (A3) 



where the Lagrange derivative of & in (|A2|) has been ex- 
pressed explicitly in (|A3[) and v from (|A1[) has been inserted 
into the last term of (|A2[) . 

The factor db/dt represents the rotation of essentially 
radial unit vector 6 with angular velocity Q and therefore 
the first term has the magnitude of Qc. The second term 
in eq. (|A3f) represents the contribution to acceleration that 
results from curvature of magnetic field lines and is zero at 
the dipole axis (solely this term fully represents the original 
model of the RC69). The third term has the same magnitude 
as the first one: \cQ x &| ~ Sic. The last term has the mag- 
nitude of Q times the local corotation velocity v CIt = X r 
and it is negligible at low altitudes where « C rt «C c. Thus, 
\a\ ~ 2Q,c and p rt ~ c 2 /a ~ R\ c /2. 

For a ^ (and C, = a) the two non-zero terms in eq. lA3l 
(and therefore a) are smaller by a factor of sin a, so that p rt 
is larger by 1/ sin a. The zero-curvature regions are located 
closer to the dipole axis by the sin a factor. However, the 



delay-radius relation retains its usual form (a- independent), 
because the sin a factor is cancelled by the 'not a great circle' 
effect. 

Yet another derivation of eq. |T} can be found in 
Thomas & Gangadgara (2005) (see their eq. 29). 



APPENDIX B: ABERRATION AND 
PROPAGATION TIME DELAYS 

The aberration is understood as the change of emission di- 
rection of radio waves caused by the transition from the 
Lorentz frame instantaneously comoving with a corotating 
emission point to the inertial observer frame (IOF). I as- 
sume that in the corotating frame (CF) the radio waves are 
emitted in the direction k' which is tangential to the direc- 
tion of dipolar magnetic field (F = ±B/\B\). Then, in the 
IOF the emission occurs in the direction: 



k = 



k' + [ 7 crt + (7crt - l)Q3crt • fcQ/ftrtl&rt 
7crt(l+/3crt -k') 



(Bl) 



where 7 crt = (1 - /?f rt ) 1/2 , and /3 crt = v CTt /c = Q x f/c. 
Eq. (IB1[) is just the general Lorentz transformation for an 
arbitrary velocity vector v, in which I substituted k — v/c 
and k — if' /a 

For low emission altitudes (r R\ c ) and near the dipole 
axis (small angle approximation) eq. (|F31|) can be consider- 
ably simplified because k! is not far from radial and v C it is 
small and purely azimuthal. Using the spherical coordinates 
(r, (f>, 9), with 5-axis along H and with the azimuth <f) increas- 
ing in the direction of the corotation velocity iT cr t = x r 
the aberration formula becomes: 



{lAk,e k ) « (l,<f>' k ,+r/R lc ,e' k , 
= k' + (0,r/i?ic,0), 



(B2) 



which means that in our specific case the aberration ba- 
sically rotates the emission direction forward in phase by 
r/R\ c rad, regardless of the emission colatitude (DRH). Note 
that eq. (|B2|) refers to the unit vectors of emission direction 
after they have been parallel-shifted to the origin of the co- 
ordinates' frame: k — (l,<j>k,9k), k' = (1, <j>' k i , 9y)- 

Propagation time delays (also called 'retardation') refer 
to the travel time of radio waves within the pulsar magneto- 
sphere: radiation emitted deeper in the magnetosphere will 
be detected later than the simultaneous emission from high 
altitudes. The delays are taken into account through the ref- 
erence to the center of the neutron star (NS). Let us consider 
a point located at the radial position r at the instant t em 
when it is emitting towards an observer. Under normal con- 
ditions such a point is located closer to the observer than 
the NS center. Radiation emitted from this point will there- 
fore be detected earlier by r(t em ) -n bs/c than an imaginary 
signal simultaneously emitted from the center of the NS. 
Here r(t em ) -n bs is a projection of r (taken at the instant of 
emission) on the observer's line of sight (defined by the unit 
vector n ba) and I assume that the radio waves travel at the 
speed of light c (dispersion is ignored) . Note that the propa- 
gation time advance/delay is not independent of aberration: 
to find r(t em ) at the moment of emission, we first need to 
know the aberrated emission direction k(r). However, since 
the angle between r(t e m) and n bs is small, in the lowest 
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order approximation the propagation time advance simply 
becomes: r(i em ) • n bs/c ~ r/c and the associated phase ad- 
vance is r/Ri c , just as in the case of aberration. 

This paper has been typeset from a Tp^X/ I^Tf^X file prepared 
by the author. 
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